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We study the dimensional crossover from 2D to ID type behavior, which takes place in the 
thermal excited rippling of a graphene honeycomb lattice, when one of the dimensions of the layer is 
reduced. Through a joint study, by Monte Carlo (MC) atomistic simulations using a quasi-harmonic 
potential and analytical calculations, we find that the normal-normal correlation function does not 
change its power law behavior in the long wavelength limit. However the system size dependency 
of the square of out of plane displacement (/i^) changes its scaling behavior when going from a 
layer to a nanoribbon. We show that a new scaling law appears which corresponds to a truly ID 
behavior and we estimate the ratio of the sample dimensions where the crossover takes place as 
R2D^iD ~ 1.609. Having explored a wide number of realistic systems sizes, we conclude that 
narrow ribbons present stronger corrugations than the square graphene sheets and we discuss the 
implications for the electronic properties of freestanding graphene systems. 
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Introduction. The very existence of graphene^, a first 
truly two dimensional crystal, is surprising because crys- 
talline order is not expected to be possible in a two di- 
mensional world. This is a consequence of a Mermin- 
Wagner theorem^ which forbids the ability of any con- 
tinuous symmetry breaking, the translational one in this 
case, in two dimensions. However, graphene exists be- 
cause the atoms can explore the other dimension by mov- 
ing perpendicular to the sheet direction. In fact, the 
stability of a two-dimensional crystalline membrane em- 
bedded in a three dimensional world has been studied ex- 
tensively in connection to biological membranes^. It has 
been concluded that the anharmonic coupling between 
the in plane and the out of plane phonons stabilizes an 
asymptotically fiat phase but strongly corrugated by the 
presence of ripples, whose manifestation is the divergence 
of the mean square of the displacement perpendicular to 
the sheet in the thermodynamical limit. The recently 
experimental demonstration of the existence of intrinsic 
ripples- in graphene has revived this subject of study. 
Moreover, numerical simulations using a realistic C-C in- 
teratomic potential supports the adequacy of the scaling 
theory of membranes in the continuum medium approach 
for graphene^^ The connection of these ripples and their 
fluctuations with the electronic transport properties is a 
subject of great interest nowadays^. 

In addition to graphene sheets, ribbons where electrons 
are confined in the nano scale in the transversal direction, 
has been the subject of a lot of interest. The reason is 
that the lateral confinement may induce the presence of 
a gap in the electronic spectrum which is essential for 
many future nanoelectronics applications. However, the 
current success for the synthesis of the semiconducting 
graphene-based nanoribbons is still far below the expec- 
tations. Lithographic patterning of graphene into nan- 



odimensions has difficulties in controlling the nanopat- 
tern size and edge qualities. It is expected that these dif- 
ficulties will be overcome in the near future. Therefore 
prediction of the structural stability as well its interplay 
with the electronic transport of graphene nanoribbons is 
strongly desirable. 




FIG. 1: Esquematic view of the corrugated graphene hon- 
eycomb lattice. The dimensions of the flat sample are given 
by Lx — SlxCiQ and Ly — \/3lyao, where ao = 1.42 A is the 
graphene lattice constant and, lx (ly) is the number of cells in 
the armchair (zigzag) -type edge. The total number of Carbon 
atoms is = 4:lxly. For this example /cc = 3 and ly — 10 and 
hence AT = 120. 

In this letter, we study the structural stability of 
graphene ribbons, extending the previous analysis con- 
ducted for the case of a square layer- . From the elasticity 
theory and using the Self Consistent Screening Approx- 
imation (SCSA) to account for the coupling between in 
and out of plane phonons, we analyze the dimensional 
crossover which takes place when one of the dimensions of 
a layer is reduced. The normal-normal correlation func- 
tion retains its power law behavior in the long wavelength 
limit. However, the mean square of the height fluctua- 
tion diverges with the length of the ribbons signaling an 
instability for the case of thin ribbons. We check this 
prediction by numerical simulations with the Montecarlo 
method using a phenomenological quasiharmonic poten- 
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tial. We conclude that ribbons should be strongly more 
rippled as compared with layer. 

Stability condition from the elasticity theory. Let 
h{x^y) be the measure of the out-of-plane displacement 
of a particle sited at the (x, y) point, in a flat configura- 
tion, of a membrane. Then, the total bending energy is 
given by 



Eq = dx I dy K (y h{x^y)) 



(1) 



where k is the bending rigidity. We are interested in 
comparing the case = Ly which represents a 2D mem- 
brane with the one Ly » corresponding to a ribbon 
running in the y-direction. It is convenient to bring eq.([T]) 
to the momentum space formulation 



(Ax, Ay) 

Eo= ^ (|k|')^|Mk)|' 

k=-(A^,A„) 



(2) 



where k = (kx-,ky) is the two dimensional momentum 
vector given as /c^ = f^^i, i = x, 7/ in terms of the integer 
summation indexes . Note that this supposes peri- 
odic boundary condition (PBC) in both directions. We 
will also use PBC in the Montecarlo simulations we will 
show in the next Section. Therefore, we take into account 
only the bulk excitations within the ribbon. The values 
of the width that make the ribbon unstable toward the 
thermally excitated rippling should be taken as an upper 
limit of stability. Edge excitations, which are not taken 
into account in our approach, could destroy even thicker 
ribbons than which we consider here. In eq. (|2]) Kx{^^y) 
is the short distance ultraviolet cutoff of the order of the 
inverse of the lattice constants ax{ciy). The prime over 
the summation sign indicates that the (0, 0) point should 
be excluded because it corresponds to a rigid transla- 
tion. The size dependency we are interested in will be 
related to this low momenta limit that appears in all the 
k-summation. 

The stability of this system can be analyzed according 
to the size dependency oi < > the mean square of the 
height and < 9'^ >^< \Vh\^ > the one of the bending 
angle formed by the normal to the surface at a point with 
the global z-axis. For the isotropic membrane {L = Lx = 
Ly) we obtain for these quantities 

1 

<h^ >=— / dkkH{k) - L^ < >- logL (3) 

47r 2tt 

where we have introduced the height-height correlation 
function given by H{k) =< /i(k)/i(— k) >= q^^- 

Note that this harmonic treatment predicts a logarith- 
mic divergence with the system size which is interpreted 
as a crumpled instability. Moreover, as we discuss lat- 
ter this instability is cured by anhamonic interactions 
with the in plane phonons. Before this, let us compare 



with the situation of a thin ribbon {Ly » Lx). In this 
case the spectrum of the bending phonons uol^^ = tvk."^ = 

is dominated by the Ux = 



term. Equation ^ now takes the form 



< >= 



4L. 
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The instability here is of course more severe that in the 
2D case. This is the first indication that thin ribbons 
are much more sensitive to thermal out plane fiuctuation 
than a square membrane. 

However, the bending modes are coupled to the in- 
plane ones. A simple way to include this coupling is by 
taking into account that in the elastic two dimensional 
theory given by the energy 



pi 



(5) 



the strain tensor Uij should include a quadratic term in 
h (A is the lame constant and ji the the shear modulus). 
Their lowest order is therefore given by 



1 



[{djUi + diUj) + dihdjh] 



(6) 



Ui are the x and y Cartesian components of the displace- 
ments from the equilibrium situation. 

The simplest but otherwise precise treatment as com- 
pared with the Montecarlo numerical simulation^ is 
the Self Consistent Screening Approximation (SCSA) of 
LeDoussal and Radzihovsky^. The in-plane phonons are 
integrated out in the partition function, generating an 
effective quartic interaction between the out of plane 
modes. The correlation function i^(k) acquires a self- 
energy correction H~^{k) = nk^ S(k) which is self- 
consistently determined. The resulting coupled integral 
equations could be analytically solved in the long wave- 
length limit under the assumption ^ ^(k) = Zk^~'^ 
Z is a non universal amplitude and 7^ 0.821 is obtained 
by solving analytically the integrals in terms of the V 
function. A crucial point for the present analysis is that 
77 cannot not depend on the system size (at least for big 
enough systems), because the integral in k involved in its 
calculation is convergent with the low momenta infrared 
limit. Therefore, we expect that 77 will not change when 
we go from a square membrane to a ribbon. This will be 
checked by Monte Carlo(MC) simulations. 

Coming back to the membrane and including the renor- 
malization of i^~^(k) given by the SCHA, Eq. (j3j) be- 
comes 



dk 



L^-"?, < 6»2 >^ L-"^ (7) 



The bending angle converges for L — > oo implying a long 
range order on the normals. However < h? > diverges 
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which is interpreted as the origin of a strong ripphng 
of the membrane which in fact observed in suspended 
graphene samples. 
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FIG. 2: (a) (Color online) q^H{q)/N against at T = 300 
K for different 2D system sizes as indicated on the plot. The 
dashed lines show the asymptotic behaviors, (b) (/i^) as a 
function of the average linear system size L = ^ LxLy com- 
pared to the scaling law (/i^) = C2dL^~'^ with C2D = 0.0026 
and rj — 0.85 (dashed line). Both axes are in logarithmic 
scale. This Figure has to be compared with the figs. 2 and 3 
of Ref . i6|- 



Monte Carlo simulations for the honeycomb lattice. 
The scaling laws described above have been obtained by 
analyzing the asymptotic behavior due to thermal fluc- 
tuations of continuous membranes. Recently, its validity 
was proven on a discrete 2D lattice, the graphene layer^. 
This was done with MC simulations using a potential 
suited for the Carbon compounds by adopting PBC on 
square samples. In fig. [2] we reproduce these results as 
a previous step before the study of the ribbons. Figure 
[2] (a) shows q^H{q)/N , which is the normal-normal cor- 
relation function because the unit vector normal to the 
surface n is connected with by n ^ — Vh. It is 
calculated over 2D layers of several sizes as indicated on 
the plot. Besides this, fig. Efb) shows {K^) as a function 
of the average linear system size L = y^L^^L^ compared 
to the eq. ([7]). These results were obtained through ordi- 
nary MC simulations including the wave-moves and usin^ 
the quasi-harmonic (QH) potential introduced in Ref. 6. 
The details of the calculation, including the implementa- 
tion of the wave- moves, follow the lines described in this 
reference. We have found the right acceptation range 
in the simulations by using amplitudes for wave-moves 
of the type A = C^jm/ L^''''''^ with C^jm = 0-01 and 
Rscai = 0.32. The observed agreement should be taken 
as a check of the reliability of our sampling procedure. 

From the previous section it is expected that when one 
of the dimensions of the square sample is reduced a dis- 
tinct behavior should appear in the scaling of {h'^)^ oth- 
erwise, q^H{q)/N should show the same scaling. Figure 
[3] (a) shows that the asymptotic behavior of q^H{q)/N 



remains valid starting from the square sample with 
N=37888 by reducing Lx as indicated (the same behav- 
ior is obtained by decreasing Ly). In fig. [3] (b) it can 
be observed that the scaling behavior of {K^) against L, 
corresponding to three different samples, is followed only 
for certain widths but it is not fulfilled when the sample 
becomes narrow enough. In all the cases we have verified 
that the same behavior is obtained by decreasing Ly. 

As previously discussed we have adopted PBC in all 
the simulations. In this form we avoid edge reconstruc- 
tion effects on the samples which could lead to unwanted 
effects for our study. We have limited the width of the 
studied ribbons to be large enough values in a way that 
bulk effects accounted by PBC are an important part of 
the physics of this problem. 
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FIG. 3: (a) (Color online) q^H{q)/N against g at T = 300 K 
calculated for different ribbon widths L^c, starting from the 
2D AT = 37888 {L^ = 315.24A)layer, as indicated on the plot. 
The asymptotic behavior remains valid independently of the 
ribbon width as was discussed previously in the text . (b) {h^) 
as a function of the average linear system size L — LxLy 
compared to the scaling law (/i^) = C2dL^~^ used in fig|2](b). 
Open symbols correspond to square lattices, being square for 
N = 37888, triangle for N = 19440 and diamond for N = 
8640. Filled symbols represent ribbon type lattices obtained 
from the square ones by reducing the Lx dimension. Starting 
from each 2D layer, in all cases the scaling law is broken at a 
certain value that we analyze further below. Both axes are in 
logarithmic scale. 

The one- dimensional behavior. From fig. [3] (a) it can 
be established that the value of 77 does not change when 
going from the square membrane to a ribbon. Now, a 
similar analysis leading to Eq. [7] can be undertaken for 
the thin ribbons. In this case, the mean square fluctua- 
tions acquire a dependence toward the long direction Ly 
and the scaling laws become 



y 



(8) 



Lx Lx 
The main result of this work is present in fig. Hlwhere it 
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FIG. 4: (Color online) {h^) versus the longitude of the ribbon 
Ly for several ribbons width Lx as indicated. The dotted line 
serves to compared to the scaling law {{h^)) — {Cid / Lx)L^~^ 
with CiD = 0.0055. The dashed line represent the 2D scaling 
for the widest ribbon (see the text). The crossover 2D to ID 
can be seen clearly in this case. Both axis are in logarithmic 
scale. 

can be observed the scaling of {K^) against Ly measured 
with the MC simulations for different ribbon widths as 
indicated. The comparison with the eq. (|8]) using the 
law {Cid/ Lx)Ly~^ with Cid = 0.0055 shows an excel- 
lent agreement in all the cases for the larger values of Ly 
where the ID aspect is stronger. The lowest Ly point of 
the Lx = 153.36 A sample (star symbol) corresponds to 
the 2D lattice N=37888 where, as we already discussed, 
the expectation value of fluctuations is described by Eq. 
([7]). For this reason we also include on the plot this law 
using the fitted value C2D = 0.0026, and Lx = 153. 36A. 
This point is the only one included in this figure which 
corresponds to a 2D layer. The change of the scaling, 
thus, emerges clearly and the fittings allow us, by equal- 
ing both laws, Eq. (j8j) (ID) and Eq. ([7j) (2D), to es- 
timate the ratio R2D^id in which the crossover should 
take place as 

T ^ 2/(77-4) 

i?2DolD = ^ = ^ 1.609 (9) 



Conclusions. The long range order on the ribbons is 
not restored by the inclusion of 77 since eq. ([8]) implies 
that < 0^ > remains unbounded in the thermodynami- 
cal limit and therefore the free thin ribbons are not sta- 
ble against a crumpled transition. For the finite sam- 
ples studied here this transition has not appeared. In 
the experimental situation, the free standing samples are 
always supported at the edges, hence the induced ten- 
sion can provide room for a higher stabilization of the 
ribbonsiSi. 



As the QH potential employed here is softer than those 
ones best suited for graphene^^^^^, to extend our conclu- 
sions to the graphene nanoribbons we emphasize that, 
R2D^iD could be greater than the estimated here. How- 
ever we do not expect a qualitative change in the overall 
behavior which is related to general scaling arguments. 
Even more, the use of more realistic Carbon-Carbon po- 
tentials would be desirable in the studies where the edge 
reconstruction is taken into account in addition to ther- 
mal excited ripples or the more stable configurations of 
low number of Carbon atoms^^. Finally, we remark that 
R2D^iD do not represent the limit in which the crumpled 
instability would disassemble the lattice, but instead, 
where the crossover takes place followed by the change 
on the scaling laws described here. A strong corrugation 
could imply big changes on the electronic spectra of the 
nanorribons and may even result in gap openings on the 
electronic spectra^ 
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